Nonlinear coupling of continuous variables at the single quantum level 
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We experimentally investigate nonlinear couplings between vibrational modes of strings of cold 
ions stored in linear ion traps. The nonlinearity is caused by the ions' Coulomb interaction and gives 
rise to a Kerr- type interaction Hamiltonian H — h\n r n s , where h r , h s are phonon number operators 
of two interacting vibrational modes. We precisely measure the resulting oscillation frequency shift 
and observe a collapse and revival of the contrast in a Ramsey experiment. Implications for ion trap 
experiments aiming at high-fidelity quantum gate operations are discussed. 
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Kerr nonlinearities in optical media caused by disper- 
sively acting third-order susceptibilities play an impor- 
tant role in nonlinear optics. There is considerable inter- 
est in realizing a Hamiltonian H = hxn r n s coupling two 
field modes by a cross-Kerr nonlinearity as this interac- 
tion might be used for quantum-nondemolition measure- 
ments of single photons [l|, for the realization of quan- 
tum gate operations in photonic quantum computation 
0, H,rJ, HI and for photonic Bell state detection Un- 
fortunately, the required strong nonlinearities are diffi- 
cult to realize experimentally 0, . 

In trapped ion quantum computing, continuous quan- 
tum variables occur in the description of the joint vi- 
brational modes of the ion string. The normal modes 
are of vital importance for all entangling quantum gates 
as they can give rise to effective spin-spin couplings in 
laser- ion interactions |9J . The normal mode picture nat- 
urally appears when the ion trap potential is modelled as 
a harmonic (pseudo-)potential and the mutual Coulomb 
interaction between the ions is linearized around the ions' 
equilibrium positions [lJJ]. In this way, the collective ion 
motion is described by a set of independent harmonic 
oscillators with characteristic normal mode frequencies. 
Small deviations from this picture arise because of the 
nonlinearity of the Coulomb interaction giving rise to a 
cross-coupling between the normal modes. While non- 
linear Coulomb couplings between ions found lots of at- 
tention in ion t rap experiments investigating order-chaos 
transitions [ll|, [ijj , its effect in the normal mode regime 
has mostly gone unnoticed. Here, we are not interested 
in resonant mode coupling [13j | but rather in dispersive 
cross-Kerr effects leading to shifts of the normal mode 
frequencies. In this paper, we show that the Coulomb 
nonlinearity gives rise to a Kerr-type Hamiltonian and 
we quantitatively measure the strength of the induced 
frequency shift. This frequency shift is of importance 
for ion trap experiments aiming at pushing the fidelity 



of entangling quantum gates towards the fault-tolerant 
threshold. 

For two cold ions of mass m and charge e held in an 
anisotropic harmonic potential characterized by trap fre- 
quencies oj z ,lo x = ujy = uij_, with uj z < u>±, the equilib- 
rium positions of the ions are given by n = (0, 0, ±zq), 
where zq = {e 2 / {lQirtomu 2 )) 1 / 3 . Separating the center- 
of-mass (COM) and the relative ion motion by introduc- 
ing the coordinates R = (r 2 + i"i)/2 and r = (r 2 — ri)/2 
yields the potential energy 

Upot = UC° M (R) + mw z 2 (z 2 + 4(* 2 + V 2 ) + 

' z ^ 

Expanding the Coulomb potential to fourth order around 
(0, 0, z ) by setting z = z + u and keeping only those 
terms giving rise to a cross-coupling between normal 
modes, we find 

U;t\ = mu 2 u 2 + mcu r 2 (x 2 +y 2 ) + V® + V™ , 

where w s = \/3uj z , uj r — J lo 2 ^ — lo\ are the stretch and 
the rocking mode frequencies [3] and where V^,V^ 
defined by 

y(3) = H^l u{x 2 + y 2^ yW = _^l u 2 {x 2 + 8j 

zo z£ 

describe the cubic and quartic nonlinearities. Next, we 
use perturbation theory for a calculation of the energy 
shifts e(n s , n%,n v r ") of the quantum states described by the 
quantum numbers (rt s ,n^,n^) where n s and n% ,y specify 
the number of stretch mode and rocking mode phonons. 
In first-order, we calculate only the shift induced by the 
quartic term as the cubic term does not contribute. In 
second order, the quartic term can be neglected. The 
shift of the stretch mode frequency by the rocking mode 
quantum numbers is given by 



2 



with the cross-Kerr coupling constant 

\ 4lu£ — bjj J \uj r J \OL A mc A ) 

where a denotes the fine structure constant and c the 
speed of light. This expansion is valid as long as Lo r /uj s 3> 
(xq/zq) and \2uj r — lu s \/lu s 3> (uq/zq) where xq (uq) 
are the spatial half width of the rocking (stretch) mode 
ground states, respectively. The resonance 2to r = co s cor- 
responds to a parametric coupling converting two rock- 
ing mode phonons into a single stretch mode excitation. 
In experiments with thermally occupied rocking modes, 
the cross-mode coupling causes dephasing of the stretch 
mode motion. 

In our experiments, two 40 Ca + ions are confined in 
a linear Paul trap with radial trap frequencies of about 
u)±/2tt — 4 MHz. By varying the trap's tip voltages 
from 500 to 2000 V, the axial center-of-mass frequency 
uj z is changed from 860 kHz to 1720 kHz. The ions are 
Doppler-cooled on the 51/2 Pi/2 transition. Side- 
band cooling on the Si/ 2 <-> D5/2 quadrupole transition 
[To] ] prepares the stretch mode in the motional ground 
state |0) s . Simultaneous cooling of stretch and rocking 
modes is accomplished by alternating the frequency of 
the cooling laser exciting the quadrupole transition be- 
tween the different red motional sidebands. Motional 
quantum states are coherently coupled by a laser pulse 
sequence exciting a single ion on the \S) = Si/2{m = 
— 1/2) \D) = D b / 2 { m = ~ 1/2) transition with a fo- 
cused laser beam on the carrier and the blue sideband. 
Internal state superpositions (\S) + e l *|D))|0) can be 
mapped to motional superpositions |-D)(|0) + e l *|l)) by a 
7r pulse on the blue motional sideband and vice versa. We 
discriminate between the quantum states S1/2 and -D 5 / 2 
by scattering light on the 5*1/2 P1/2 dipole transition 
and detecting the presence or absence of resonance flu- 
orescence of the individual ions with a CCD-camera. A 
more detailed account of the experimental setup is given 
inRef. [IE[l3. 

For a measurement of the stretch mode coherence, 
a Ramsey experiment between motional states is per- 
formed by a laser interacting with only one of the ions, 
the second ion being just a spectator that modifies the 
normal mode structure. Starting from the state |5)|0), 
the superposition state |-D)(|0) + |1)) is created by a n/2 
pulse on the carrier transition followed by a tt pulse on the 
blue sideband. During a waiting time of duration t, the 
state evolves into |-D)(|0) 4- e*^|l)), where <j> is a random 
variable. Finally, the motional superposition is mapped 
back to superposition of internal states (\S) + e l ^\D))\0) 
and probed by a tt/2 pulse on the carrier transition fol- 
lowed by a quantum state measurement. We measure the 
coherence by varying the phase of the last tt/2 pulse from 
to 2tt and measuring the contrast C (r) of the resulting 
Ramsey pattern as a function of r. The measurement 
shown in Fig. Q] exhibits a strong reduction in contrast 
after only 2 ms. This fast decay could be attributed 
neither to motional heating (because the stretch mode 
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FIG. 1: Contrast C(r) of a Ramsey experiment on the stretch 
mode as a function of the Ramsey waiting time r. The axial 
trap frequency was u> z = (27r)1486 kHz, the radial modes 
were cooled close to the Doppler limit. The initial contrast 
is limited by imperfect sideband cooling, spurious excitation 
of the second ion by the focussed laser beam and magnetic 
field noise. The subsequent loss of contrast is caused by a 
dephasing of the stretch mode oscillation. The inset shows 
the pulse sequence of the Ramsey experiment with C denoting 
carrier pulses and B blue sideband pulses. 
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FIG. 2: Contrast C(r) of a spin echo experiment on the 
stretch mode as a function of the Ramsey waiting time r. 
Insertion of the spin echo increased the coherence time by 
nearly two orders of magnitude as compared with a simple 
Ramsey experiment. 



heating rate was less than 0.02 phonons/ms), nor to an 
instability of the trapping potential (since a coherence 
time of about 50 ms was measured on the center-of-mass 
mode). To probe the dynamics of the dephasing mech- 
anism, a spin echo experiment was performed on a mo- 
tional superposition |0) + For this, the populations 
of the |0) and |1) states were exchanged in the middle of 
the experiment by a carrier 7r-pulse sandwiched between 
two blue sideband 7r-pulses. As can be seen in Fig. [21 
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FIG. 3: (a) Shift A(j> of the spin echo phase pattern by a 
rocking mode phonon observed by measuring the D state 
population of ion 1 as function of the phase of the last spin 
echo pulse. The symbol (o) labels experiments with an extra 
phonon in the second spin echo time, (o) experiments with no 
extra phonon. (b) Single phonon shift of the spin echo phase 
pattern as a function of the waiting time. The frequency shift 
caused by a single rocking mode phonon is inferred from the 
slope of the function A(j>(r). 



now it takes the spin echo contrast C(r) about 100 ms 
to decay to 50% of its initial value - a clear indication 
that the stretch mode frequency is fairly stable over the 
duration of a single experiment but randomly changing 
from experiment to experiment. This behaviour is ex- 
pected for a dephasing caused by a thermally distributed 
rocking mode phonon number that takes on random val- 
ues at the start of each experiment. To confirm this de- 
phasing mechanism, we remeasured the Ramsey contrast 
C(r) with both rocking modes cooled close to the ground 
state and observed an increase of the coherence time by 
a factor of 20 as compared to the experiment shown in 
Fig. [T] 

Having established the cross-coupling between the 
stretch mode and the rocking mode as the dephasing 
mechanism, we are interested in quantifying the cross- 
coupling strength by measuring the stretch mode fre- 
quency shift caused by a single rocking mode phonon. 
In principle, such a measurement could be performed by 
preparing the rocking modes in a n r = or n r — 1 Fock 
state and subsequently measuring the stretch mode oscil- 
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FIG. 4: Shift of the stretch mode frequency pattern by a single 
rocking mode phonon as a function of the axial trap frequency 
uj z . The solid line is the frequency shift predicted by eq. ijl). 
the dashed line shows only the quartic contribution, and the 
dotted line is a fit to the data using a heuristic scaling law 
dv a /dn r oc cjf with (5 — 3.25. 



lation frequency in a Ramsey experiment. This, however, 
would require ultra-stable high-voltage sources for keep- 
ing the axial trap frequency ui z stable to within 10~ 6 c<j 2 
or less. To make the experiment more robust against 
technical noise, we carried out a spin echo experiment 
instead, where the other ion was excited on the blue 
sideband of one of the rocking modes directly before the 
start of the second spin echo waiting time r (see inset 
of Fig. [3ja)). If this pulse excites the second ion into 
the -D5/2 state, there is exactly one additional phonon 
in the rocking mode during the second half of the spin 
echo whose frequency shift is not compensated by the 
echo sequence. We adjust the duration of this blue side- 
band pulse so that the excitation is successful in 50% 
of the experiments and sort the spin echo data into two 
classes according to the quantum state of the second ion 
at the end of the experiment. By measuring the phase 
shift A0(t) between the two data sets, we are able to 
infer the shift caused by a single rocking mode phonon. 
This procedure does not even require cooling the rocking 
modes to the ground state. Fig. (3^a) shows the resulting 
phase shift for a waiting time r = 30 ms. For a calcula- 
tion of the single phonon frequency shift dis s /dn r , we plot 
the phase shift as a function of the spin echo time r, fit a 
straight line to the data and determine its slope. For the 
data measured at lu z /2it = 1716 kHz shown in Fig.[3Jb), 
the frequency shift is dv s /dn r = 20.5 Hz/phonon. This 
measurement procedure was carried out for different axial 
center-of-mass voltages lo z while keeping the transverse 
oscillation frequency lo±_ fixed. In Fig. U the frequency 
shift is plotted as a function of lu z . The experimentally 
measured frequency shift is slightly smaller than the shift 
predicted by perturbation theory. The disagreement be- 
tween experimental data and theoretical model is cur- 
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FIG. 5: Collapse and revival of the contrast in a Ramsey experiment probing the motional coherence of the two lowest stretch 
mode quantum states. The revival occurs at the time predicted by the spin echo experiment for a frequency lu z /2tv — 1716 kHz. 



rently not understood. It cannot be attributed to the 
dynamical nature of the trap potential. Using a refined 
pseudopotential [l8| leads only to marginal corrections. 

For a thermally occupied rocking mode, the change 
induced in uj s is always an integer multiple of the single 
phonon shift. Therefore, after the initial collapse of the 
Ramsey contrast (Fig. Q} a revival [l9[ is to be expected 
for a time r* = 2ir/x = 2ir(duj s / '<in r ) _1 . Operating the 
trap at uj z /2t: — 1716 kHz, we indeed observed a revival 
at the predicted time r*. In this experiment, one of the 
rocking modes was cooled to the ground state while the 
other one was prepared in a thermal state with n r « 9 
(lcr-confidence interval {5,17}). Fig. [5] shows a fit to the 
data using the function 

C(t) = e - 7r |(e l *™" T )| = e-~< T \n r + 1 - n r e lXT \-\ 

assuming a thermally occupied rocking mode and an 
overall loss of contrast cx e _7r accounting for techni- 
cal noise and motional heating. A fit to the data yields 
a revival time r* = 50.5(5) ms, an average vibrational 
quantum number n r = 9 (2) that is consistent with 
an independent measurement of n r and a decay rate 
7 = 0.004(3)s _1 . The existence of revivals is a fur- 
ther proof of the quantized nature of the rocking motion. 
The small loss of contrast shows that the probability of 
a change in the rocking mode phonon number within the 
interval [0, t*] is quite low. 



The observed frequency shifts need to be taken into 
account in quantum gate realizations operating on the 
stretch mode. For the parameters [2(| used in ref. 0, 
eq. (QJ predicts shifts as big as 100 Hz/phonon, giving 
rise to a loss of fidelity of about 0.1% for n = 1. 

In summary, we have investigated a nonlinear quantum 
effect giving rise to a cross-coupling of harmonic oscil- 
lators that can be described by a Kerr-like Hamiltonian 
H (x n s n r . In a two-ion crystal, the nonlinearities lead to 
a dephasing of the relative ion motion that manifests it- 
self as a collapse of the Ramsey contrast that revives once 
all oscillations get in phase again. While the nonlinearity 
is fairly small for the trap parameters investigated here, 
it could be made bigger by tuning the normal mode fre- 
quencies closer to the resonance u) r = 2oj s so that it might 
be used for creating entangled motional states of these 
oscillators. On the other hand, for quantum computing 
experiments aiming at fault-tolerant quantum gates, the 
observed nonlinearity points to the necessity of cooling 
all spectator modes to the ground state or working with 
transversally very stiff ion traps. 
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